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Abstract 


The  basic  concepts  of  geometrical  optics  together  with  the  additional  assump¬ 
tions  that  lead  to  the  "geometrical  optics  approximation"  are  described  here.  The 
eikon&l  equation  is  derived  and  the  relationship  of  exact  electromagnetic  theory  in 
the  limit  as  X  —  0  to  geometrical  optics  is  made  evident.  The  application  of  the 
"geometrical  optics  approximation"  to  phase  analysis  and  synthesis  is  described 
and  an  example  of  synthesis  is  presented.  The  concept  of  power  flow  in  ray  tubes 
is  used  to  obtain  approximations  to  power  distributions  in  the  antenna  aperture,  in 
the  focal  region,  and  in  the  far  field.  Ray  analysis  is  used  to  determine  those  feed 
locations  in  the  focal  region  that  will  most  nearly  collimate  the  far-field  rays  that 
lie  in  certain  desirable  planes.  The  Theorem  of  Malus  is  used  to  formulate  the 
equal  path  length  law  and  applications  are  given.  Focal  surfaces  (or  caustics) 
relative  to  a  rectilinear  congruence  are  defined  and  then  used  to  present  a  geometrical 
optics  description  of  the  focal  region.  The  equations  of  the  focal  surfaces  of  a  para¬ 
boloid  receiving  a  plane  wave  20°  off-axis  are  calculated  and  photographs  of  three- 
dimensional  models  of  the  focal  surfaces  are  shown. 
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Application  of  Geometrical  Optics  to  the  Design 
and  Analysis  of  Microwave  Antennas 

1.  INTRODUCTION 

This  report  is  concerned  with  the  laws,  principles,  and  procedures  of  geomet¬ 
rical  optics  that  are  applicable  to  the  design  and  analysis  of  microwave  antennas. 
Geometrical  optics,  considered  as  a  zero  wavelength  approximation  to  exact  electro¬ 
magnetic  wave  theory,  is  very  accurate  in  the  design  ard  analysis  of  optical  focusing 
devices  because  optical  wavelengths  are  extremely  small  compared  to  the  aperture 
dimensions  of  optical  systems.  At  microwave  frequencies,  the  wavelength  is  not 
always  relatively  small  compared  to  the  aperature  dimensions  of  microwave  sys¬ 
tems.  Geometrical  optics,  however,  although  certainly  an  approximation  is  still 
sufficiently  accurate  to  produce  meaningful  and  useful  results,  even  lor  antennas 
with  aperture  dimensions  as  small  as  five  wavelengths.  The  advantages  of  certain 
microwave  components  over  existing  optical  components  and  the  relaxation  in 
mechanical  tolerance  requirements  due  to  the  finite  wavelength  allow  the  exploitation 
of  geometrical  optics  analysis  and  design  at  microwave  frequencies  in  certain  cases 
considerably  beyond  that  achievable  in  optics.  For  example,  low-loss  isotropic 
artificial  dielectrics  in  a  wide  range  of  index  of  refraction,  phase  and  amplitude 
control  of  sources  and  receivers,  and  aspheric  as  well  as  nonrotational  symmetric 
reflecting  and  refracting  surfaces  are  all  available  to  the  antenna  designer. 

(Received  for  publication  5  May  1967) 
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In  genera!,  geometrical  optics  is  concerned  with  the  analysis  and  synthesis 
of  optical  systems  to  the  approximation  that  diffraction  and  interference  can  be 
neglected.  In  an  isotropic  medium,  classical  geometrical  optics  assumes  that  the 
power  flows  along  paths  called  rays  at  a  velocity  characteristic  of  the  medium  and 
that  there  exists  a  family  of  surfaces  known  as  wavefronts  that  are  everywhere 
normal  to  the  rays.  Point-to-point  correlation  between  wavefronts  can  be  estab¬ 
lished  by  the  rays  and  no  power  is  assumed  to  be  present  in  regions  where  there 
are  no  rays.  It  is  evident,  that  if  all  the  wavefronts  are  given,  all  the  rays  are 
determined  and  vice  versa.  Classical  geometrical  optics,  therefore,  neglects  wave¬ 
length,  phase,  and  the  vector  nature  of  electromagnetic  wave  motion.  For  micro- 
wave  applications,  it  is  most  useful  to  extend  the  classical  theory  to  include  effects 
of  the  above -neglected  factors.  The  form  of  the  extension  is  justified  by  the  asymp¬ 
totic  solution  as  w  — -  oo(X  —  0)  of  the  exact  electromagnetic  field  equations  to  be 
discussed  in  Section  3.  The  extension  consists  of  introducing  wavelength  as  a  small 
but  finite  quantity,  identifying  the  wavefronts  with  equiphase  surfaces,  and  at  each 
point  on  a  ray  in  a  homogeneous  medium  introducing  the  electromagnetic  field  vec¬ 
tors  _E  and  H  and  relating  them  as  in  a  plane  wave  propagating  along  the  ray.  We 
shall  refer  to  this  extension  as  the  geometrical  optics  approximation. 


2.  THE  EiKONAL  AND  THE  EIKONAL  EQUATION 

Using  the  geometrical  optics  approximation,  assume  for  a  particular  wave 
motion  in  a  scurce-tree  isotropic  medium  that  the  equiphase  wavefronts  are  given 
by  the  level  surfaces  of  the  function 

L  -  Ux,  y,  z) 

and  that  the  phase  6  on  a  general  wavefront  W  is  given  by 

4>  =  wt  -  -  L(x,  y,  z)  (1) 

where  u  is  the  angular  frequency,  c  is  the  velocity  in  free  space,  and  (x,  y,  z) 
is  any  point  on  W.  The  function  L(x,  y,  z)  is  known  as  the  eikonal  and,  together 
with  the  wavefront  velocity  in  the  medium,  it  completely  describes  the  given  wave 
motion  from  the  standpoint  of  classical  geometrical  optics. 
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Consider  a  general  ray  path  C  given  by  the  equations 

x  =  x(s) 

y  =  y(s) 

z  =  z(s)  (2) 

where  s  is  arc  length  measured  along  C  from  a  fixed  reference  wavefront  WQ  , 

given  by  L(x,  y,z)  =  0  to  a  propagating  wavefront  W  (see  Figure  1).  Let  the 
position  of  W  along  C  as  a  function  of  time  be  given  by  the  formula 

s  =  s(t)  (3) 


where  s  increases  with  .  The  wavefront  propagation  velocity  v  along  C  is 
then 


(4) 


The  value  of  L(x,  y,z)  on  the  wavefront  W  will  change  as  W  propagates  along  C. 
Since  time  and  position  along  C  are  related  by  (2)  and  (3),  L(x,  /,  z)  may  be  con¬ 
sidered  as  a  function  of  either  s  or  t  .  Thus 


|vl| 


dL 

dL  _  dt_ 
ds  ~  ds 
dt 


(5) 


Since  the  phase  4>  on  the  propagating  wavefront  W  is  constant,  differentiating 


Figure  1.  Ray  Pat-is  and  Wavefronts 
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where  n  is  the  index  of  refraction  of  the  medium.  The  partial  differential  equa¬ 
tion  (7)  satisfied  by  the  eikonal  L(x,  y,  z)  is  known  as  the  eikonal  equation.  * 

The  eikonal  equation  can  be  used  to  determine  the  curvature  of  ray  paths  in  an 
inhomogeneous  medium.  Assume  a  wave  motion  in  an  inhomogeneous  medium  pre¬ 
scribed  by  the  eikonal  L(x,  y,  z)  and  let  C  be  a  ray  path  with  arc  length  s  in  the 
direction  of  the  wave  motion.  At  each  point  on  C  let  S  be  the  unit  tangent  vect&r 
in  the  direction  of  s  ,  let  flpbe  the  unit  principal  normal  vector,  and  let  p  be  the 
radius  of  curvature.  The  vector  VL  evaluated  at  any  point  P  is  normal  to  the 
wavefront  passing  through  the  point  P.  Hence,  at  each  point  on  C  ,  the  vector 
VL  is  along  C  and  in  view  oi  me  eikonal  equation  (7), 


VL 

n 


(8) 


The  first  Frenet  formuk  (Hildebrand,  1949)  of  differential  geometry  can  be  written 


fi.  =  (§  •  V)  g  =  -  g  X  (VX  g)  =  4?  •  (9) 

ds  “  —  —  -  P 

Substituting  (8)  into  (9)  we  have 

fjl  =  -gX  |v(±)xvl]  =  g  X  [ V (in  n)Xg] 

and 


1_ 

P 


flp  •  V(ln  n) 


-  fi 
n  ~p 


Vn  . 


(10) 


Equation  (10)  shows  that  the  ray-path  curvature  ^  is  related  to  the  rate  of 
change  of  the  index  of  refraction  n  notmai  to  the  ray  path.  In  particular,  for  a 
homogeneous  medium  (n  =  constant)  the  curvature  is  zero  and  the  ray  paihs  are 
straight  lines. 


i 

3.  GEOMETRICAL  OPTICS  AS  A  ZERO  WAVELENGTH  APPROXIMATION 

For  a  source -free  nonconducting  isotropic  homogeneous  medium.  Maxwell's 
equations  for  the  electric  field  E  and  the  magnetic  field  H  are 
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VXE  =  -jw/iH  (11a) 

VXH  =  jco«E  (11b) 

V.  E  =  0  (11c) 

V-  H  =  0  (lid) 


where  w  is  the  angular  frequency,  t  is  the  permittivity  of  the  medium,  p  is  the 
permeability  of  the  medium,  and  an  time  dependence  is  assumed.  Equations 
(11a,  b,  c)  combine  to  produce 

V2  E+  k2  E  =0  (12) 


where  k  =  is  the  phase  constant  of  the  medium. 

To  obtain  a  zero  wavelength  or  high  frequency  approximation,  the  electric 
field  is  assumed  in  the  form  of  an  asymptotic  series  (Luneberg,  1944;  Kline  and 
Kay,  1965;  Kouyoumjian,  1965)  in  descending  powers  of  w  as  follows: 


E  (x,  y,  z,oj) 


e-jk0L(x,  y,  z) 


£ 


(jtf)”1 


(13) 


where  E^x,  y,  z)  is  real  and  kQ  =  <->V  «0P0  is  the  phase  constant  of  free  space. 

Substituting  (13)  into  (12)  and  independently  equating  to  zero  che  coefficients  of 
2 

w  and  to,  the  highest  powers  of  u  present,  yields 

|vl|2  =  n2  (14) 


and 

2. 

(VL-V)E  +X-^-E  =  0,  (15) 

where  n  =  ^ —  is  the  index  of  refraction  of  the  medium. 

Note  tha*?  (14)  is  the  eikonal  equation  and  hence  L(x,  y,  z)  is  the  eikonal.  As 
to  becomes  large,  the  leading  term  of  (13)  predominates  and  becomes  E^p  »  the 
high  frequency  approximation  to  the  electric  field,  as  follows: 

-jk  L(x,.y,  z) 

Ehf  =  S^fx,  y,  z)e  .  (16) 

Since  E  (x,  y,  z)  is  real,  the  equiohase  surfaces  L(x,  y,  z)  =  constant  of  E„_ 
are  identical  with  the  wavefronts  of  geometrical  optics.  Again,  the  rays  can  be 
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defined  as  the  family  of  curves  normal  to  the  wavefronts.  In  the  present  case  of 
an  isotropic  homogeneous  medium,  the  rays  are  straight  lines  (see  Section  2). 

At  each  point  on  a  ray  9  is  a  unit  vector  tangent  to  the  ray  [see  Eq.  (8)  ] . 

Thus, 

VL-V  =  n^p-V  =  nS  -  V  =  n^  (17) 

where  is  the  directional  derivative  in  the  direction  of  a  ray  with  respect  to 
arc  length  s  along  the  ray.  Equation  (15)  can  now  be  written 


_1  V_1 

+  2  n 


0  . 


(18) 


The  solution  to  ( 18)  along  a  ray  through  the  point  ( xQ,  yQ,  zq)  can  be  expressed  in 
the  form 


E0<s>  =  E0(s0)  e 


s 

■hi 


v2l 


ds 


(19) 


where  so  is  the  distance  from  a  reference  wavefront  to  the  point  ( xQ,  yQ,  zQ).  Thus, 
Eq(s)  at  any  point  on  a  ray  is  completely  determined  once  its  value  is  known  at  one 
point  on  the  ray.  This  is  a  very  clear  statement  of  the  geometrical  optics  property 
of  point-to-point  correlation  from  wavefront  to  wavefront  along  a  ray.  It  is  also 
evident  from  (19)  that  the  direction  of  Eq  is  the  same  for  all  points  on  a  ray  (ex¬ 
cept  oossibly  for  sense). 

2 

Substituting  (13)  into  (11c)  and  equating  to  zero  the  coefficient  of  u>  ,  the 
highest  power  of  w  ,  leads  to  the  result 


VL  •  Erf  =  0  (20) 

and  hence  at  each  point  on  a  ray  is  no.  n.al  the  ray.  Substituting  (13) 

into  (11a)  and  retaining  only  the  highest  order  terms  in  u>  yields  H^..  i‘*c  high 
frequency  approximation  to  the  magnetic  field,  as  follows: 

Huc<  =  -—2-  VLXE  e"jkoL  r-l/L’s  X  EUI,  .  (21) 

-HF  -jcjw  -o  \J  v  —  —  HF 

Equations  (20)  and  (21)  show  that  EHF  and  along  a  ray  are  related  precisely 

as  in  a  plane  wave  propagating  in  the  direction  of  the  ray. 


1 .  .'tv Z  J :>?&*£- 


Simplification  of  (19)  is  best  accomplished  bv  considering  the  vector  field 


F  =  Kn  3  ,  K  = 


where  at  each  point  R.  and  R_  are  the  principal  normal  radii  of  curvature  of  the 

z  j 

wavefront  passing  through  that  point  (see  Appendix  A).  The  quantity  K  =  ^  ^  - 
evaluated  at  a  point  on  a  surface  is  known  as  the  Gaussian  curvature  of  the  surface 
at  that  point.  It  can  be  shown  that  the  surface  integral 


ffr  ■  «  .  0 


where  T.  is  any  closed  surface  lying  in  the  isotropic  homogeneous  medium  being 
considered.  The  divergence  theorem  then  requires  that 

V  .  F  =  0 

throughout  the  medium  and  hence 

(n  3  •  V)  K  --  -K  V  •  n  3  =  -K  V2L  . 


in  view  of  (17)  it  follows  that 


f  • 


Hence, 


m’  .  e 
K7-TJ 


/4 


Equation  (19)  can  now  be  written 


SoM  -  wv-^Sh j  • 

o 


By  applying  (17)  to  the  eikonal  L(x,  y,  z)  and  using  (14)  we  have 


2  dL 
=  n  =  n  -t— 
ds 


VL  •  VL 


$$$£  .2"-’ i\ ^&r  K  k4-  AJLT5?*#  w*Mf  tr***- 
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or 


dL  =  nds . 


For  a  homogeneous  medium  n  is  constant  and 


L  =  ns  +  LQ  . 


(?5) 


Substituting  (24)  and  (25)  into  (16)  yields 


E  -  £  (s  )J  —te)  e"jk°Lo  e"jks 

-HF  -o'soW  K(s  )  e  e 

o 


(26) 


Note  that  the  amplitude  factor*/  is  precisely  equivalent  to  the  geometrical 

IT(A\  ™  rw'sr»/ 


K(d)  1  *m-0. 

optics  power  factor  ]^g)  shown  in  (B3)  of  Appendix  B. 


4.  FERMAT’S  PRINCIPLE  AND  SNELL’S  LA»S 

The  optical  path  length  (OPL)  along  a  ray  patn  C  is  defined  as  th°  line  inte¬ 
gral  Jc  nds  where  n  is  the  index  of  refraction  of  the  medium  and  s  is  arc  length 
along  C.  Fermat's  principle  states  that  electromagnetic  energy  traveling  between 
two  points  will  follow  any  ray  path  that  makes  the  OPL  integral  stationary.  Clearly 
then,  ray  paths  in  a  homogeneous  medium  will  be  straight  lines. 

Snell's  laws  of  reflection  and  refraction  of  rays  at  a  boundary  surface  between 
two  different  media  can  be  derived  directly  from  Fermat's  principle,  but  will 
merely  be  stated  here.  The  law  of  reflection  states  that  the  incident  ray  and  the 
reflected  ray  lie  on  the  same  side  of  the  boundary  surface,  are  coplanar  with  the 
normal  to  the  boundary  surface  at  the  point  of  reflection,  and  make  equal  angles 
i  with  the  normal.  Thus,  0.  =  where  0.  is  the  angle  of  incidence  and  0f  is 
the  angle  of  reflection,  as  shown  in  Figure  2.  This  reflection  law  can  be  expressed 
in  any  of  the  equivalent  vector  forms 


n  X 

<£r-£i> 

=  o  , 

(27a) 

n  • 

(-r  +  *i> 

=  o  , 

(27b) 

5  = 

— r 

«0  f 

CM 

1 

Cfl  { 

•  n)  fi  , 

(?7c) 

5.  = 

— l 

5  -  2(5 

—r  — r 

*  n)  n  , 

(2'7d) 

9 


where 


and  S.  is  in  the  direction  of  the  indicent  ray,  $r  is  in  the  direction  of  the  reflected 
ray,  and  n  is  normal  to  the  boundary  surface  (see  Figure  2). 


'TRANSMITTED  RAY 


INDEX  OF  REFRACTION  n, 


Figure  2.  Reflection  and  Refraction  ai  a  Boundary 

The  law  of  refraction  states  that  the  incident  ray  and  the  transmitted  or  refracted 
ray  lie  on  opposite  sides  of  the  boundary  surface,  are  coplanar  with  the  normal  to 
the  surface  at  the  point  of  refraction,  and  satisfy  the  condition 

n.  sin  0 .  -  n,  sin  0 .  ,  ,  h  I  (28) 

l  l  t  t  ( 

where  n.  and  n{  are  the  indices  of  refraction  in  the  incident  and  transmitted 
regions  respectively  and  6  .  and  0 1  are  the  angles  of  incidence  and  refraction 

respectively  (see  Figure  2).  This  refraction  law  can  be  expressed  in  the  vector 

form 

n  X(n.  S.  -  n  §  )  =  0  (29) 

—  1  ~ *1  I  -H 

where  n  ,  S.  ,  n.  ,  and  n,  are  as  defined  above  and  S,  is  a  unit  vector  in  the 

—  —l  i  t  -t 

direction  of  the  transmitted  ray. 

Combining  (27a)  and  (29),  we  have  as  a  statement  of  both  of  Snell’s  laws 

n  X  n.  S.  =  fiXn  S  =  n  X  n,  5,  (30) 

—  l  —l  —  r  — r  t  — t 

where  nf  =  n..  .  that  is,  the  incident  and  reflected  regions  have  the  same  index  of 
refraction. 
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5.  PHASE  ANALYSIS  AND  PHASE  SYNTHESIS 


I 


The  geometrical  optics  approximation  can  be  applied  usefully  to  certain 
problems  in  phase  analysis  and  phase  synthesis.  An  application  to  phase  synthesis 
is  given  in  this  section. 

If  L(x,y,  z)  is  the  eikonal  function,  then,  according  to  (1),  the  phase  4>  at  a 
point  P  with  coordinates  (x,  y,  z)  is 


<t>  =  u>t  -  “  L(  x,  y,  z)  . 


If  Wq  is  the  wavefront  L(x,  y,  z)  =  0  then  the  phase  <pQ  over  the  whole  wavefront 

W  is 
o 


I 

f 

i 

i 


=  ut  • 

The  phase  at  P  relative  to  the  pnase  on  W 
between  P  and  W  ,  will  be 


o 


*  =  <?  - 


L(x,  y,  z) 


that  is,  the  phase  difference  <p 

(31) 


Assuming  that  the  ray  through  P  intersects  WQ 
to  Section  4,  the  optical  path  length  from  P  to 


OPL 


nds  , 


at  point  Pq  ,  then,  according 
P  is  defined  a4 

0  / 

/ 


where  n  is  the  index  of  refraction  of  the  medium  and  the  path  of  integration  is 
along  the  ray  from  PQ  to  P  .  The  phase  difference  <t> ^  between  P  and  Pq 
S  can  be  expressed  in  terms  of  OPL  as  follows: 


where  is  the  free-space  wavelength.  Noting  that  ^  =  ^pi-and  compar¬ 
ing  (31)  with  (32),  it  is  evident  that  the  eikonal  function  L(x,  y,  z)  evaluated  at  P 
is  equal  to  the  OPL  from  the  reference  wavefront  to  P  .  The  geometrical  optics 
determination  of  phase  will  clearly  fail  in  any  region  where  two  or  more  rays  pass 
through  each  point. 

Now  consider  the  problem  of  designing  a  point-source-fed  reflector  in  two  dimen¬ 
sions  in  a  medium  of  index  n  that  will  synthesize  to  within  the  geometrical  optics 

/ 

/ 

/ 
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approximation  a  given  phase  distribution  along  a  line  (Sletten,  1958).  Let  the  point 
source  be  located  at  the  point  F  with  coordinates  (0,  a)  and  let  the  x  axis  be  the 
line  along  which  it  is  desired  to  achieve  the  phase  function  f  (x)  .  If  a  ray  reflected 
from  the  reflector  at  the  poiut  P  with  coordinates  (x,  y)  intersects  the  x  axis  at 
the  point  PQ  with  coordinates  (xq  ,  0),  as  shown  in  Figure  3,  then  it  follows  from 

(32)  that 

(-^r)n  Vx2  +  (y-a)“  +(2f)  n  ^(x-xq)2  +  y2  =  -f(xo)  +  Cj  (33) 

where  C  ^  is  a  constant  that  can  be  adjusted  to  control  the  position  of  the  reflector 
and  the  extent  of  useable  aperture.  For  xq  fixed  (33)  as  a  function  of  x  and  y 
represents  an  ellipse  with  foci  at  F  and  Pq  .  With  xq  variable  (33)  represents 
a  one  parameter  family  of  ellipses  that,  in  general,  has  an  envelope.  To  each 
ellipse  in  the  family  there  corresponds  a  unique  point  on  the  envelope  where  the 
ellipse  is  a  tangent  to  the  envelope.  By  using  the  envelope  as  the  reflector  surface, 
it  is  clear  that  at  each  point  on  the  envelope  the  incident  ray  is  precisely  reflected, 
as  from  the  ellipse  corresponding  to  that  point,  and  hence  arrives  at  the  x  axis  in 
proper  phase. 

The  envelope  of  the  family  of  ellipses  is  given  by  the  simultaneous  solution  of 

(33)  with  the  partial  derivative  of  (33)  with  respect  to  xq  ,  that  is, 

=  f'(xo)  .  (34) 


A 


A 

2wn 


Figure  3.  Geometry  of  Reflector  to 
Produce  Desired  Phase  Along  a  Line 
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and 


B  *  ^  »'<«„>  • 

An  experimental  line  source  feed  designed  by  the  above  procedure 
for  off-focus  feeding  of  a  paraboloid,  was  very  successful  in  producing  a  well- 
focused  off-axis  fan  beam. 


6.  POWER  FLOW  IN  RAY  TUBES 

One  of  the  important  assumptions  of  geometrical  optics  is  that  the  power  flows 
along  the  ray  paths.  Therefore,  a  ray  path  diagram  presents  an  overall  picture  of 
the  power  flow  for  a  lossless  source-free  medium  to  within  the  accuracy  of  this 
approximation.  The  totality  of  rays  that  pass  through  any  given  closed  curve 
constitutes  a  ray  tube;  and,  under  steady  state  conditions,  the  total  power  flowing 
across  any  cross  section  of  a  ray  tube  must  be  constant  since  no  power  can  flow 
across  the  lateral  surfaces  of  the  cube.  Thus,  as  a  ray  tube  cross  section  decreases 
the  power  density  increases  and,  conversely,  as  a  ray  tube  cross  section  increases 
the  power  density  decreases.  These  two  cases  correspond  to  converging  and  diver¬ 
ging  rays  respectively.  In  particular,  if  dA  is  a  differential  area  on  a  wavefront 
W  'and  if  the  ray  tube  passing  through  dA  intersects  the  wavefront  W*  in  the 
differential  area  dA1  (see  Figure  4),  then,  the  total  power  flow  through  dA  must 
equal  the  total  power  flow  through  dA  1  .  Thus,  we  have 

P  dA  =  P1  dA  ’  (36) 


Figure  4.  A  Ray  Tube  Between 
Wavefronts 


where  P  is  the  power  density  at  dA 
and  P1  is  the  power  density  at  dA’  . 

It  is  clear  that  the  ray  tube  concept  of 
power  flow  will  break  down  at  focal 
points  since  at  such  points  the  ray  tube 
cross  section  vanishes  and  (36)  predicts 
infinite  power  density. 

The  concept  of  power  flowing  in 
ray  tubes  as  applied  to  the  region  be¬ 
tween  the  primary-feed  and  the  aper¬ 
ture  is  particularly  useful  in  the  design 
and  analysis  of  antennas.  In  far -field 
considerations,  it  is  also  useful  in  the 
design  and/or  analysis  of  certain  shaped 


» »  *  *•***>  & ' 
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beam  and  off-focus-fed  antenna  systems.  The  ray  tube  concept  finds  little  appli¬ 
cation  in  far-field  considerations  relative  to  well-focused  pencil  beam  systems 
since  in  such  cases  the  pattern  is  determined  entirely  by  diffraction. 


7.  POWER  DISTRIBUTION  IN  THE  APERATURE 


In  using  geometrical  optics  to  investigate  antenna  aperture -plane  power  dis¬ 
tributions  or  to  design  primary  feeds  for  optimum  illumination,  the  important 
quantity  furnished  by  the  ray  tube  concept  is  an  approximation  to  the  relative 
power  distribution. 

To  determine  the  relative  power  distribution,  it  is  necessary  to  know  the 
absolute  power  distribution  to  within  only  an  arbitrary  multiplicative  constant. 
Hence,  if  either  the  absolute  or  the  relative  power  distribution  is  known  on  a  wave- 
front  W  (see  Figure  4),  then,  in  using  (36)  in  the  form 


to  determine  the  geometrical  optics  approximation  to  the  relative  power  distribu- 

dA 

tion  on  W'  »  it  is  necessary  to  know  ^7  to  within  only  a  multiplicative  constant. 

For  a  well-focused, point-source-fed  antenna  system,  the  usual  procedure  for 
determining  the  geometrical  optics  approximation  to  the  relative  power  distribution 
on  the  antenna  apert-  re  plane  for  a  given  primary  feed  power  pattern  is  to  first 
analyze  che  system  when  fed  by  an  isotropic  point  source.  Ail  wavefronts  will  be 
spheres  for  the  isotropic  feed  and  the  power  density  will  be  constant  over  any 
fixed  wavefront  W  .  Also,  any  differential  element  of  area  dA  on  the  wavefront 
W  will  be  proportional  to  the  differential  solid  angle  dft  subtended  by  dA  as 
measured  at  the  feed.  A  general  ray  trace  of  the  antenna  system  will  determine 
the  aperture  coordinates  of  the  exit  ray,  as  functions  of  the  ray  direction  coordinates, 
as  measured  at  the  feed.  These  relations  will  determine  the  area  dA'  -'n  the  aper¬ 
ture  plane  W1  (assumed  to  coincide  with  a  wavefront)  that  corresponds  to  the  dif¬ 
ferential  element  of  area  dA  in  the  wavefront  W  .  Since  dft  =  KdA,  where  K 
is  a  constant,  the  quantity  is  proportional  to  -^r  and  can  therefore  be  used  in 

the  calculation  of  the  relative  power  distribution  on  W'.  If  the  primary  feed  has  a 
relative  power  pattern  P (0,<p),  where  0  and  <p  are  the  ray  direction  coordinates, 
then  the  geometrical  op*ics  approximation  to  the  relative  power  distribution  on  the 
aperture  plane  is  proportional  to  P (6 ,  <p)  ,  where  (in  theory)  6  and  <t>  can 

be  expressed  in  terms  of  the  aperture  coordinates. 

The  geometrical  optics  approximation  to  the  relative  power  distribution  on  the 
aperture  of  a  wcil-focused.line-source-fed  cylindrical  antenna  system  can  be 
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determined  by  procedures  analogous  to  those  used  for  the  point-source  case.  It 
will  be  assumed  that  the  relative  power  pattern  of  the  line-source  feed  is  of  the 
form  P (0 ,  z)  =  g  j (0  )©2  ’  w^ere  z  1S  linear  coordinate  parallel  to  the  feed 

line  and  0  is  the  angular  coordinate  around  the  feed.  Since  the  antenna  system 
is  cylindrical,  all  rays  from  any  one  point  on  the  feed  line  remain  in  a  plane  nor¬ 
mal  to  the  feed  line  and,  in  all  such  planes,  ray  trace  diagrams  are  identical. 

Thus,  power  will  only  flow  parallel  to  these  planes  and  the  identity  of  the  ray 
trace  diagrams  guarantees  that  the  relative  power  distribution  on  the  aperture 
plane  in  the  z  direction  will  be  proportional  to  gg (z)  .  This  reduces  the  problem 
by  one  dimension  and  leads  to  the  consideration  of  the  power  flow  in  a  typical  plane 
normal  to  the  feed  line.  For  a  feed  isotropic  in  0  ,  the  wavefronts  are  circles  and 
the  power  density  on  these  circles  is  uniform.  If  ds  is  a  differential  element  of 
arc  length  on  a  fixed  wavefront  C  and  if  d0  is  the  differential  angle  at  the  feed 
subtended  by  ds,  then,  since  C.  is  a  circle,  it  follows  that  us  is  proportional 
to  dfl  .  A  gensral  ray  trace  will  relate  0  to  h  ,  where  h  is  the  exit  ray  aperture 
coordinate  normal  to  z  ;  this  relation  will  determine  the  differential  element  ah 
in  the  aperture  plane  (assumed  coincident  with  a  wavefront)  that  corresponds  to  ds  . 
Since  ds  -  K'  d0  where  K'  is  a  constant,  it  follows  that  is  proportional  to 
-gp-  and  can  therefore  be  used  in  the  calculation  of  the  relative  power  distribution 
on  the  aperture  plane  with  respect  to  the  h  coordinate.  For  the  given  feed  with 
relative  power  pattern  P(0,  z)  =  gj(0)g2(z)  >  the  geometrical  optics  approxima¬ 
tion  to  the  relative  power  distribution  on  the  aperture  plane  is  proportional  to 

jrt 

g1(0)g2(z)  ,  where  (in  theory)  0  can  be  expressed  in  terms  of  h. 

As  a  simple  example  of  this  procedure  applied  to  a  cylindrical  system,  con¬ 
sider  a  parabolic  cylinder  reflector  of  focal  length  f  with  its  focal  line  coincident 


y 


! 
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/ 
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Figure  5.  Central  Plane  Section  of  a 
Parabolic  Cylinder  Antenna 


with  the  z  axis.  Figure  5  shows  a 
central  plane  section  of  this  reflector 
with  polar  coordinates  (r,  0)  describ¬ 
ing  the  reflector  surface  and  the  linear 
coordinate  h  designating  position  in 
the  aperture  plane  x  =  0  .  The  equa¬ 
tion  of  the  reflector  surface  in  (r,  0 ) 
coordinates  is 


1  +  cos  0 


f  sec  (0  /  2). 


The  primary  feed  located  at  F  (see 
Figure  5)  is  assumed  to  radiate  uni¬ 
formly  in  0  and  all  rays  originating 
at  F  are  assumed  to  be  reflected 
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parallel  to  and  in  the  direction  of  the  negative  x  axis.  Thus, 
h  =  r  sin  0  =  2ftan(6/2) 
dh  =  f  sec2  (8  /2)dS  rd0  ■ 

and 


dd  1 
dh  r  ‘ 

If  the  given  primary  feed  has  a  power  pattern  P(0,  z)  =  g,(#)g,,(z),  then,  P  , 

1  Ct  £L 

the  geometrical  optics  approximation  to  the  relative  power  distribution  on  the 
aperture  normalized  at  the  center  of  the  aperture,  is  given  by 


j_  g1(g)g2(2) 

Pa  "  r  g1<0)g2(0) 


(38) 


This  is  certainly  a  correct  result  from  the  geometrical  optics  standpoint  since  it 
clearly  demonstrates  that  the  power  spreads  radially,  that  is,  inversely  propor¬ 
tional  to  distance,  only  in  planes  normal  to  the  feed  line  and  only  over  the  dis¬ 
tance  from  the  feed  to  the  reflector,  that  is,  over  the  distance  r  .  After  reflec¬ 
tion,  the  rays  are  collimated  and  no  more  spreading  takes  place  out  to  the  aperture 
plane. 


8.  POWER  DISTRIBUTION!  IN  THE  FOCAL  REGION 

In  the  analysis  of  an  imperfectly^focused  antenna  suen  as  a  spherical  reflector 
or  an  off-axis-fed  paraboloid,  it  is  frequently  desired  to  obtain  an  estimate  of 
the  relative  power  distribution  along  some  focal  surface  or  curve  in  the  focal  re¬ 
gion.  In  this  case,  the  ratio  of  a  differential  element  of  surface  area  of  the  incoming 
plane  wavefront  to  the  corresponding  differential  element  of  surface  area  or  arc 
length  in  the  focal  region  is  an  estimate  of  the  relative  power  distribution  on  the 
appropriate  focal  surface  or  curve. 

Consider,  for  example,  the  problem  of  determining  the  power  distribution 
along  the  axis  of  a  spherical  reflector  of  radius  "a"  illuminated  by  a  plane  wave 
incident  along  the  axial  direction  {see  Figure  6).  By  symn  etry  all  ray::  after 
reflection  pass  through  the  reflector  axis.  The  height  h  of  an  incident  ray  and 
the  coordinate  z  of  the  intersection  of  the  corresponding  reflected  ray  with  the 
reflector  axis  are  related  to  the  angle  6  as  follows; 
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h  =  a  sin  6 


Hence, 


2  cos  6 


SPHERICAL  V 
REFLECTOR 


REFLECTOR  AXIS' 


-INCOMING  PLANE  WAVE 


REFECTOR  CENTER  OF 
CURVATURE 


Figure  6.  Central  Section  of  a  Spherical  Reflector 


.  2  2 
n  =  a 


(-£)■ 


The  differential  area  dA  generated  by  re  /olving  the  line  segment  dh  about  the 
reflector  axis  is 

dA  =  2  it  h  dh  . 

All  incident  rays  passing  through  dA  will,  after  reflection,  pass  through  the 
differential  line  element  dz.  Thus,  the  power  distribution  P  along  the  reflector 


■>  ■  ■  2,hBfr 


where  B  is  a  constant.  From  (39)  we  have 
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From  (39  )  we  have 

lim  z  =  a/2  .  (41) 

h-*  0 


This  is  known  as  the  paraxial  focus  point.  Normalizing  the  power  distribution  to 
unity  at  the  paraxial  focus  point  the  relative  power  distribution  along  the  re¬ 
flector  becomes 

3 

r-,  _  a _ 


9.  POWER  DISTRIBUTION  IN  THE  FAR  FIELD 

In  the  use  of  the  ray  tube  concept  to  approximate  antenna  far-field  power  ^ 

patterns,  the  underlying  assumption  is  that  the  radiated  power  after  it  leaves  the  f 
final  surface  of  the  antenna  flows  along  ray  tubes- not  only  up  to  the  aperture  plan& 
but  indeed  on  out  into  the  antenna  far  field.  This  assumption  is  clearly  invalid  for 
well-focused  antenna  systems,  that  is,  for  systems  all  of  whose  exit  rays  are  nearly 
parallel,  since  in  such  cases  the  pattern  is  determined  entirely  by  diffraction. 
However,  for  systems  that  are  not  well-focussed  in  one  or  in  both  planes,  the  ray 
tube  concept  provides  a  useful  approximation  to  the  relative  far-field  pattern.  In  any 
plane  in  which  the  exit  rays  are  not  parallel,  the  power  How  per  unit  directional  angle  or 
per  unit  solid  angle  in  the  exit  ray  system  is  determined  by  the  ray-tube  concept 
as  a  function  of  angular  direction;  this  function,  after  appropriate  normalization, 
is  the  geometrical  optics  approximation  to  the  far-field  antenna  pattern. 

As  an  example  of  the  use  of  this  far-field  approximation,  consider  the  design 
of  a  hne-source  fed  cylindrical  reflector  to  produce  a  fan  beam  of  prescribed 
shape  (Spencer,  1943).  Let  S  be  the  plane  of  symmetry  of  this  antenna  normal 
to  the  line  source.  The  cylindrical  geometiy  guarantees  that  all  exit  rays  lying 
in  planes  normal  to  S  are  parallel  and  therefore  S  will  be  the  plane  of  the  fan 
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beam.  It  is  desired  to  determine  the  shape  of  the  reflector  cioss  section  that  will 
cause  ray  dispersion  in  the  plane  of  S  appropriate  to  croduce  the  prescribed 
fan  beam  shape.  It  is  assumed  that  the  primary  feed  pattern  Pf(0)  and  the  de¬ 
sired  secondary  power  pattern  P(  <b)  are  given.  The  angles  0  and  <p  are  mea¬ 
sured  positive  clockwise  in  the  plane  of  S  .  (See  Figure  7a.)  For  practical  pur¬ 
poses,  the  useful  portion  of  the  primary  feed  pattern  P^'O)  lies  between  its  -10  dB 
points  and  will  be  designated  by  the  interval  0  0s0  2  •  Generally,  the  function 

P (<p)  vanishes  everywhere  except  within  some  interval,  say  <p  ,  where  it 
takes  on  a  prescribed  form.  Thus,  rays  from  the  line-feed  incident  on  the  reflec¬ 
tor  with  inclination  angle  0  ,  0  j-  (!  ^0^  must  produce  reflected  rays  with  in¬ 
clination  angle  <>,  <t>^  <0  5:02  (see  Figure  7a>. 

The  problem  now  is  to  determine  a  one-to-one  correspondence  oe  tween  the 
incident  and  reflected  rays— that  is,  a  function  <$>  -  p{0)  such  that  *he  reflected 
ray  system  has  a  power  distribution  corresponding  to  P(4>).  In  a  differential  sec¬ 
tor  d0  of  the  incident  ray  system,  the  power  flow  is  proportional  to  Pf(0)d0; 
in  a  corresponding  differential  sector  d <t>  of  the  reflected  ray  system,  the  power 
flow  is  proportional  to  P ($)  d<*>  .  Thu3, 


Pf(0)d0  =  KP  (<*>)  dip 


wht  re  K  is  determined  by  the  relation 


/ 


e. 


pf(6)  de  =  k 


P(<P)  d<t> . 


which  is  equivalent  to  requiring  equal  total  powers  in  the  incident  and  reflected  ray 
systems.  We  also  require  the  power  in  the  incident  ray  system  in  the  interval  0  j 
to  0  to  be  equal  to  the  power  in  the  reflected  ray  system  in  the  interval  <t> j  to  <t> . 
Thus, 


(42) 


and  this  expression  determines  the  required  relationship  <p  =  <p(9)  . 

We  now  turn  to  the  determination  of  the  equation  for  the  reflector  curve.  With 
the  reflector  curve  denoted  by  r  =  r(0)  we  have  from  Figure  7b 
V. 
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Thus, 

r  <•>■  *  "‘(^T4)  ■»  '  |-2-^J  d« 

or 

toK)’/°-i1^l“i dB-  «« 

vhere 

rQ  =  r(0)  . 

Equation  (43)  determines  the  desired  reflector  curve  to  within  a  radial  scale  factor 

r  .  Overall  physical  dimensions,  as  well  as  the  requirement  that  the  reflector  be 
o 

in  the  far  field  cf  the  feed,  establish  a  range  of  permissible  values  for  r  . 

Generally,  the  function  P^(  6 )  is  given  graphically  while  P(<p)  may  be  given 
analytically,  graphically,  or  numerically.  Therefore,  in  most  cases  numerical 
integration  will  be  required  to  determine  4>  -  <p(  6 )  from  (42)  and  the  resultant 
relationship  will  be  expressed  in  numerical  form.  Numerical  integration  will  be 
required  to  solve  (43)  for  the  reflector  curve  r  =  r(0). 


10.  REFLECTION  FROM  A  CONDUCTING  SURFACE 


In  an  isotropic  homogeneous  medium,  the  geometrical  optics  approximation 
assumes  that  the  energy  propagates  along  the  rays  as  in  a  plane  wave.  To  within 
this  approximation,  therefore,  the  field  vectors  E  and  H  along  a  ray  are  related 
as  follows: 


a 

S  =  EXH*=^  |e{2  S 


(44a) 

(44b) 


6  *  E  -£*H=E-H  =  0 


(44c) 


where  S  is  the  Poynting  vector  and  3  is  a  unit  vector  in  the  direction  of  the  ray 
(see  Figure  8).  The  magnitude  of  E  (and  hence  of  H  and  S)  is  prescribed  by 
(24). 
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With  the  above  assumptions,  con¬ 
sider  reflection  from  a  smooth, 
perfectly-conducting  surface  V  .  Let 
E  .  be  the  electric  field  incident  on  V 
along  the  ray  direction  ,  let  _E  be 
the  electric  field  reflected  from  V 
along  the  corresponding  ray  direction 

3  ,  and  let  fi  be  the  unit  surface  nor- 

— r  — 

mal  to  V  (see  Figure  9). 


Figure  8.  Electric  and  Magnetic  Fields 
Along  a  Ray  Path 


Figure  9.  Incident  and  Reflected  Electromagnetic  Fields 
and  Ray  Tubes 


Consider  a  differential  area  dA  on  V  .  Let  dAi  be  the  differential  cross  section 
of  the  ray  tube  incident  on  dA  and  let  dArbe  the  differential  cross  section  of  the 
ray  tube  reflected  from  dA.  Evaluating  dAj  and  dAr  at  V  ,  we  have  the  fol¬ 
lowing  relations: 


dA 


dAi  dAr 

«  •  Sif  12  *  'SJ 


Since  Snell  law,  (27b)  gives 


» 
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we  have 


and  hence  at  V 

|S||  ■  I— r| '  <«> 

The  continuity  of  the  tangential  electric  field  at  the  surface  V  requires 

nX(E.  +  Er)  =  0  (46) 

and  the  plane  wave  assumption  requires 

S.  •  E.  =  9r  •  Er  =  0.  (47) 

It  follows  from  (45)  and  (46)  that  at  V 

fl  •  (E .  -  Er)  =  0  (48) 


and 


Er  -  2  <a  •  Ej)  ft  -  Ej 


(49a) 


or 


E  =  (A  *  E.)  A  -  A X (E .  X  A)  .  (49b) 

—  r  —  —  1  —  —  —  1  — 

The  magnetic  fields  Hj  and  Hr  at  V  are  given  by  (44a)  and  ( 49a)  or  (49b). 

The  above  results  show  that, under  the  assumptions  made,the  reflection  process 
locally  at  each  point  of  V  is  equivalent  to  plane-wave  reflection  from  an  infinite 
conducting  plane  tangent  to  V  at  the  point  (Silver,  1949).  The  assumptions  of 
plane-wave  propagation  along  rays  and  local  reflection  behavior  as  from  an  infinite 
plane  lead  to  useful  approximations  of  the  polarization’ characteristics  of  reflector- 
type  antennas. 


11.  RAY  C0LLIMAT10N  AND  OFF-FOCUS .FEEDING 

Consider  a  focusing  system  fed  by  a  point  source  not  necessarily  located  at  a 
focus  point  and  let  §r  be  a  unit  vector  in  the  direction  of  the  general  exit  ray. 
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The  components  of  §r  will  be  the  direction  cosines  of  the  exit  rays,  that  is, 

§r  =  cos  a  5t  +  cos  p  $  +  cos  7  2. 

A  different  set  of  exit  rays  and  associated  direction  cosines  is  generated  for  each 
different  feed  position.  If  the  direction  cosines  are  considered  as  functions  of  co¬ 
ordinates  u  and  v  that  designate  the  rays— for  instance,  exit  aperture  coordinates— 
then  total  collimation  of  the  rays  is  characterized  by  the  conditions 


(50a) 


3  cos  a. 
3  v 


3  cos  , 
3  v 


_  3  COS  7 
3  v 


(50b) 


All  exit  rays  are  totally  collimated  and  the  above  conditions  hold  for  all  u  and  v 
in  the  on-focus  case  for  a  perfectly  focusing  system.  Overall  total  collimation  does 
not  occur  in  the  off-focus  case;  but,  there  may  be  conditions  under  which  partial  coi- 
limation  is  achieved  for  certain  families  of  rays.  Rays  will  be  considered  to  be 
pa.  tially  collimated  if  they  are  all  parallel  to  the  same  plane;  rays  that  are  all 
parallel  to  two  nonparallel  planes  will  be  totally  collimated.  Requiring  any  one 
of  the  direction  cosines  to  vanish  will  determine  conditions  for  partial  collima¬ 
tion  relative  to  one  of  the  coordinate  planes;  for  instance,  setting  cos  a  =  0  will 
determine  for  each  different  feed  positionthe  rays  that  are  parallel  to  the  plane  a  =  x/2, 
that  is,  the  yz  plane.  The  location  of  these  partially-collimated  rays  is  parti¬ 
cularly  important  in  antenna  design,  as  it  determines  those  portions  of  the  focus¬ 
ing  system  that  are  the  most  effective  in  achieving  partial  collimation. 

In  practice,  most  focusing  systems  have  a  plane  of  symmetry  U  and  all  rays 
leaving  the  feed  in  U  will  remain  in  U  .  Therefore,  these  rays  are  partially 
collimated.  If  U  is  taken  to  be  the  yz  coordinate  plane,  as  shown  in  Figure  10, 
each  ray  in  U  can  be  designated  by  the  angle  0  at  which  it  leaves  the  feed  point 

(°  .  yQ  •  *0>- 

Applying  the  geometrical  optics  concepts  of  power  flow  in  the  far  field  as 
described  in  Section  9  to  the  two-dimensional  case  of  the  rays  in  U  ,  the  power 
density  P  per  unit  angle  <p  at  a  range  R  is  given  by  the  expression 


P 


K  G(0 ) 


R 


dj> 

d0 


(51) 
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Figure  10.  Section  of  Focusing 
System  in  the  Plane  of  Symmetry  U 


where  K  is  a  constant  and  G (0)  is  the 
feed  power  pattern  ij»  U.  This  expres¬ 
sion  gives  a  rough  estimate  of  the  far- 
field  power  pattern  in  the  plane  of  U 
provided  /  0  .  If  the  exit  ray  at 
angle  <f>  -  <p  corresponds  to  the  ray 
leaving  the  feed  at  angle  9  =  0Q  ,  the 
condition  . 


j 

de  |e=eo 


=  o 


(52) 


determines  the  feed  position  that  will 
produce  local  collimation  of  the  rays  in 


D  that  are  in  the  neighborhood  of  the  exit  ray  in  the  direction  Under  these  conditions, 

geometrical  optics  predicts  an  infinite  singularity  in  the  far-field  pattern,  which  may 
be  interpreted  as  a  power  pattern  maximum  in  the  direction  <£  =  ;  but.  relative 

magnitudes  obviously  cannot  be  determined  without  the  use  of  diffraction. 


12.  CONGRUENCES  AND  THE  EQUAL  PATH  LENGTH  LAW 

A  two-parameter  family  of  curves  constitutes  what  is  known  in  differential 
geometry  as  a  congruence.  The  member  curves  of  a  congruence  are  called  genera¬ 
tors.  If  these  generators  are  straight  lines,  the  congruence  is  denoted  as  recti¬ 
linear.  Furthermore,  if  there  exists  a  surface  that  is  normal  to  all  the  generators, 
the  congruence  is  designated  as  normal.  Thus,  in  a  homogeneous  isotropic  medium, 
the  normals  to  a  wavefront— that  is,  the  rays  themselves— constitute  a  normal 
rectilinear  congruence.  It  can  also  be  shown  that  a  normal  rectilinear  congruence 
possesses  not  just  one  normal  surface  but  a  whole  family  of  normal  surfaces 
(Eisenhart,  1960,  p.  393);  these  surfaces,  of  course,  correspond  to  wavefronts 
of  the  corresponding  ray  congruence. 

Rectilinear  congruences  are  not  necessarily  normal  and,  therefore,  the  Theorem 
of  Malus  (Eis.enhart,  1960,  p.  403)  is  of  fundamental  importance  in  ray-path  analysis. 
This  theorem  states  that  if  a  family  of  rays,  initially  a  normal  rectilinear  con¬ 
gruence,  is  reflected  or  refracted  any  number  of  times  by  successive  homogeneous 
isotropic  media,  the  rays  will  continue  to  constitute  a  normal  rectilinear  congruence. 
Therefore,  for  antenna  systems  satisfying  the  conditions  of  this  theorem,  both 
incident  and  exit  wavefronts  are  guaranteed  to  exist.  SiJice  wavefronts  are  equi- 
phase  surfaces  ,  it  follows  that  between  any  two  specific1  wavefronts  the  optical  path 
lengths  along  all  ray  paths  must  be  equal;  this  is  a  mo^  important  result  and  forms 
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the  basis  of  the  equal  path  length  law  and  the  associated  design  procedures  appli¬ 
cable  to  reflection  and  refraction. 

Consider  the  problem  of  designing  a  reflector  R  passing  through  a  given  point 
PQ  that  will  reflect  a  given  incident  wavefront  W  into  a  desired  reflected  wave- 
f  ont  W  '  (see  Figure  11).  The  desired  reflector  surface  is  given  by  the  locus  of 
points  P  determined  by  the  condition 

d  +  d'  =  do  +  d^  .  (53) 

It  is  very  important  to  note  here  that  Snell's 
law  of  reflection  will  automatically  be  satis¬ 
fied  at  the  reflector  surface.  In  the  corres¬ 
ponding  refraction  problem,  it  is  desired  to 
determine  the  boundary  surface  S  between 
two  different  homogeneous  media  passing 
through  a  given  point  PQ  such  that  a  given  Figure  ^  Wavefronts  and  Rays  for 

incident  wavefront  W  is  converted  into  a  Reflection  From  a  Curved  Surface 

desired  refracted  wavefront  W  1  (see 

Figure  12).  In  this  case,  the  surface  S  is  the  locus  of  points  P  determined  by 
the  condition 

nd  +  n'd  1  =  nd  +  n'd  '  ,  (54) 

o  o 

and  Snell's  law  of  refraction  will  automatically  be  satisfied  at  the  interface. 


INDEX  OF  INDEX  OF 

REFRACTION  n  REFRACTION  n* 


Figure  12.  Wavefronts  and  Rays  for  Refraction 
Through  a  Curved  Interface 

In  both  of  the  above-mentioned  procedures,  it  is  necessary  that  the  reflecting 
or  refracting  surface  be  located  in  a  region  where  through  each  point  there  passes 
not  more  than  one  ray  of  the  incident  ray  congruence  and  not  more  than  one  ray  of 
the  reflected  or  refracted  ray  congruence.  If  this  restriction  is  violated,  it  will 


i 
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not  be  possible  to  design  a  reflecting  or  refracting  surface  that  satisfies  the  equal 
path  length  requirement  along  all  rays. 

As  an  example  of  the  application  of  the  equal  path  length  lav/,  consider  the 
design  of  a  reflector  to  correct  for  the  inherent  spherical  aberration  of  a  spherical 
reflector  (Holt  and  Bouche,  1964).  Let  the  radius  oi  curvature  of  the  spherical 
reflector  be  ta"and  let  the  center  of  curvature  be  the  origin.  With  the  x  axis  as 
the  reflector  axis,  an  incoming  plane  wave  along  the  x  axis  is  assumed;  the  yz 
plane  is  the  reference  wavefront.  The  spherical  reflector,  the  corrector,  and  the 
associated  ray  system  are  all  radially  symmetric  about  the  reflector  axis;  there¬ 
fore,  all  sections  in  planes  through  the  x  axis  are  the  same  (see  Figure  13). 


y 


Figure  13.  Spherical  Reflector  Section  in  the  xy  Plane 

To  fix  the  solution,  let  V  with  coordinates  (a,  0)  be  the  vertex  of  the  reflec¬ 
tor,  let  V  '  with  coordinates  (a',  0),  a  1  s  a/2,  be  the  vertex  of  the  corrector,  and 
let  F  with  coordinates  (b,  0),  b  ^a',  be  the  final  focus  point.  The  condition 
a's  a/2  guarantees  that  the  corrector  will  lie  in  a  region  where  only  one  ray  of 
the  incident  ray  congruence  and  only  one  ray  of  the  reflected  ray  congruence  will 
pass  through  each  point. 

Applying  the  equal  path  length  law  to  the  general  ray  QPP'F  and  the  axial 
ray  OVY’F  ,  we  have 

QP  +  Pp  +  P*F  -  OV  +  +  V'F 


or 


a(cos  9)  +  ^(x*  -  x)^  +  (y '  -  y)^  +  V(x'-  b)^  +  y’2 


2a  +  b  -  2a  ' 
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where 

x  =  a  cos  0  ,  y  =  a  sin  0  ,  and  y'  -  y  =  (x'  -  x)  tan  20  . 

Solving  for  x'  and  y'  ,  we  obtain- 

x'  =  a  cos  0  1  cos2  9  '  4a(a~a'>  cos  6  -  <a2  +  b2)  +(2a ±  b  -  2a')2]  (2  cos 20  -  1) 

4(b  cos^9  -  a  cos  9  -  a'  +  a) 

(55a) 

y.  =  +  (2x'  cos  9  -  ajfBfaO)  (55b) 

(2  cos20  -  1) 

'T’hese  are  the  coordinates  of  the  corrector  surface  in  terms  of  the  parameter  0  . 
Provided  a'  and  b  lie  in  the  indicated  ranges,  the  answer  given  above  is  valid 
and  unique.  (The  +  sign  for  the  y'  value  merely  indicates  symmetry  with  respect 
to  the  x  axis. ) 

13.  FOCAL  SURFACES 

Let  us  consider  a  wavefront  W  and  its  associated  ray  congruence.  Let  P  be 
a  point  on  W  and  let  PN  be  a  line  segment  normal  to  W  (see  Figure  14).  If  C 

K 1 

is  a  line  of  curvature  of  W  with  p  =  Rj  at  P  (see  Appendix  A),  then  in  the  near 

vicinity  of  P  the  normals  to  W  along  Cr  appear  to  converge  to  a  point  Q,  on 

iV  1  -* 

PN  on  the  concave  side  of  at  a  distance  R^  from  W.  Similarly,  if 

is  a  line  of  curvature  with  p  =  R,,  at  P  ,  then  in  the  near  vicinity  of  P  the  normals 
to  W  along  CR2  appear  to  converge  to  a  point  Q2  on  PN  or  the  concave  side  of 
Cr2  at  a  distance  R2  from  W  .  Thus,  to  each  point  P  on  the  wavefront  W 
there  is  associated  a  unique  pair  of  points  Q1  and  Q2  on  the  normal  line  PN. 

These  points  are  called  focal  points;  in  general,  the  totality  of  all  these  points 
constitutes  two  surfaces  which  are  called  focal  surfaces  or  caustics.  These  surfaces 
may  degenerate  into  such  forms  as  a  single  surface,  a  surface  and  a  curve,  a  single 
curve,  or  a  single  point.  It  is  important  to  note  here  that  the  focal  surfaces  are  the 
same  for  all  wavefronts  of  a  particular  ray  congruence. 

The  focal  surfaces  can  also  be  described  in  terms  of  the  envelopes  of  certain 
families  of  rays.  A  family  in  this  case  consists  of  the  totality  of  rays  passing 
through  a  particular  line  of  curvature  of  W  ,  _  ruled  surface  so  constituted  is 

known  as  a  principal  surface.  Two  typical  families  are  shown  in  Figure  14.  It  is 
characteristic  of  each  family  that  it  envelopes  a  curve  in  space  (Sj  and  S2  in 
Figure  14).  It  can  be  shown  that  each  of  these  space  curves  lies  on  a  focal  surface; 


28 


thus,  the  rays  of  a  congruence  are  tangent  to  the  focal  surfaces  along  these  envelope 
curves.  The  totality  of  all  the  envelope  curves  constitutes  the  focal  surfaces. 

In  any  specific  problem  regarding 


I 

I 

» 


Figure  14.  Ray  Paths  in  the  Prin¬ 
cipal  Planes 


the  determination  of  focal  surfaces, 
either  a  ray  congruence  or  a  wavefront 
will  generally  be  given.  Since  all  wave- 
fronts  can  be  determined  from  the  ray 
congruence  and  vice-versa,  these  two 
forms  of  tne  statement  of  the  problem 
are  entirely  equivalent.  In  actually 
carrying  out  the  computations  to  obtain 
the  focal  surfaces,  however,  the  form 
i  of  the  statement  of  the  problem  and  the 
functions  involved  will  determine  the 
best  solution  procedure.  Two  proce¬ 
dures  will  be  given  here— one  particularly 
applicable  when  the  wavefront  is  given 
and  the  other  applicable  when  the  ray 
congruence  is  given. 

Consider  the  problem  of  determining 
the  focal  surfaces  given  a  wavefront  W 
in  the  form 

P  =  x(u,  v)  X  +  y(  u,  v)  +  z(u,  v)  z  , 

where  P  is  a  vector  from  the  origin  to 
a  point  P  on  W  and  u  and  v  are 
curvilinear  coordinates  on  W  .  The  nor¬ 
mal  vector  to  W  at  P  is  given  by 


N  =  P  X  P  , 

_  — u  —v  ’ 

where  the  subscripts  denote  partial  differentiation  and  the  expression  is  evaluated 
at  P  .  To  each  of  the  two  lines  of  curvature  and  C^2  through  P  there 

corresponds  a  principal  unit  normal  vector.  These  unit  vectors  will  be  in  the 
direction  of  N  or  -N  ,  whichever  points  to  the  concave  side  of  the  correspond¬ 
ing  line  of  curvature.  Let  nj  be  the  principal  unit  normal  for  and  let  m, 

be  the  principal  unit  normal  for  The  principal  normal  radii  of  curvature 

Rj  and  R2  ji  W  at  P  are  given  by  the  solutions  to  the  quadratic  equation 
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(eg  -  f2)p  2  -  (Eg  -  2Ff  +  Ge)p  +  (EG  -  F2)  =  0  , 
where 


and  D2  =  EG  a  F2 


e  =  <ZUU  *  N)/D  , 
f  =  ■  W/D  , 

g  =  (Pw  •  N)/D  . 
=  N  •  N  . 


The  subscripts  in  tr..,  above  expressions  denote  partial  differentiation.  The  focal 
points  of  W  at  P  will  then  be  given  by  the  formulas 

Q,  =  P+R^j 


and 


&  -  ~  R2  2a  • 

As  P  ranges  over  W,  the  above  formulas  determine  the  focal  surfaces. 

The  problem  of  determining  the  focal  surfaces  for  a  given  ray  congruence  is 
now  considered.  A  ray  congruence  is,  in  general,  specified  by  giving  the  ray 
direction  at  each  point  on  a  i  eference  surface— the  reference  surface  not  necessarily 
being  a  wavefront.  Let  the  reference  surface  V  be  given  in  the  form 

R  =  >:(u,  v)  5t  +  y{u,  v)  £  +  z(u,  v)  2  ,  (56) 

where  R  is  a  vector  from  the  origin  to  a  point  R  on  V  .  Let  the  ray  directions 
at  each  point  on  V  be  specified  by  the  unit  vector  function 

w  =  w  (u,  v)  . 

If  L  is  the  ray  through  R  ,  the  distances  ctj  and  <*2  from  R  along  L  to  the 
focal  surfaces  are  given  by  the  solutions  to  the  quadratic  equation 

(EG  -  F2)  /32  +  (Eg.-  2Ff  +  Ga?  {S  +  (eg,-  f  s  f  0  s  i57) 


^  :'v‘; " 
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<c,*  -  r,-  v\\  * 


'  A  «  i.  ■-  ,'  '•  ■*  ^ 


§  "a  xto-  y  +  a  $  , 
~r  x  —  y  ~  z  — ‘ 


where 


-  2r  siny  sin  0cos  <S-  4rf 


r2  +  4f2 


2  2 
2r  sin  7  sin  0-  4rf 


r  +  4f“ 


-  sin  7 


2  2 

„  4rf  sin  7  sin  0-  8f  cos  y  , 
a  =  - 1 — .  r  - -  +  cos  y  . 

7  +  4i 


4^2  _2  ~ 

R  =  r  cos  4  8  +  r  sin  A  8-t  -7- —  §.  and  w  -  S  , 

—  _  i  4  f  —  „^r* 


the  equations  (60),  (61),  and  (62)  constitute  a  normal  congruence  in  the  usual  ;orm 
of  a  reference  surface  v/ith  ray  directions  specified  uL  each  point  on  the  reference 
surface  as  c.  function  of  curvilinear  coordinates.  The  second  procedure  201  deter¬ 
mining  the  focal  surfaces  is,  therefore,  the  most  applicable  m  this  case. 


The  fundamental  quantities  become 


9  9  2  9 

E  =  lef  (cos *7  4-  sin  y  sm  0}/A“ 


90  2  -2 

=  (16rf  sin  7  sin  0  +  8r  f  sir.  7  cos  7  cos  0)/A 


G  =  4r2  r2  +  4f“  -  (2f  sin  7  sin  0  -V  r  cos  y)2  /A2 


j £.  =  2(r  sin  7  sin  0  -  2f  cos  y)/A 


f  =  0 


g  =  2r~(r  sin  7  sin  0  -  2f  cos  7)/ A 


where 


=  r^f2 


< 
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Substituting  into  <f>7)  Had  solving  f^r  p  (with  much  labor),  we  obtain 


ft  =  — - —  +-4f2 -  |  [(r2  -  4^)  sin2y  +  8f 2 

iol  (2f  ccs  y  -  r  sin  y  sin  rf>)  ( 

-  4rf  sin  y  cos  y  sin  p  ±  (r2  +  4f2)  sin  yl 


/•}  2  To 

,  (r  -  4f  )  cos  y  +  4rf  sin  y  sin  <p 

f  t  - 


Equations  (59a  end  b)  then  give  the  focal  surfaces  in  terms  of  the  curvilinear  coor¬ 
dinates  (r,  4>). 

In  terms  of  cartesian  coordinates  (x,  y,  a)  on  the  paraboloid  surface,  the  equa¬ 


tions  for  the  focal  surfaces  become 

£  =  { x-ff}  * +  jy -®  [J 


(2f  -  z)  sin  y 


l  J  +  0X11  t  \  o 

.  i  2?  2f  cf/i  y  -  y  sin  y  j  (  *- 


+  |  z  -  B  j  1  -  - 1  (  2 

J  I  2f  cos  y  -  y  sin  -y  J  1  — 


where 


Y  i .  c°s  if~yzsinyj 


B  =  2f  -  z  +  cos  y  (z  cos  y  -  y  sin  y)  ±  (2f  -  z)  sin  y 


Calculations  for  the  particular  case  y  =  20®  have  been  carried  out  and  the 
focal  surfaces  calculated  and  constructed  (Parke  Mathematical  Laboratories, 
Report  No.  1,  1952).  Figures  16a,  b,  and  c  show  the  focal  surfaces  individually 
and  superimposed. 

14.  GENERAL  COMMENTS  ON  FOCAL  SURFACES 


The  normals  to  a  wavefront  through  a  line  of  curvature  tend,  in  general,  to 
focus  on  one  or  the  other  of  the  two  focal  surfaces.  The  physical  extent  of  the 
focal  surfaces  is,  therefore,  a  measure  of  the  focusing  ability  of  the  ray  congruence 
associated  with  the  wavefront.  If  the  focal  surfaces  are  extensive  and  far-  apart, 
the  focusing  is  poor;  if  the  focal  surfaces  are  small  and  closely  spaced,  the  focus¬ 
ing  is  good.  If  the  focal  surfaces  degenerate  to  a  line  or  to  a  point,  the  focusing 
may  be  considered  as  perfect,  depending  on  the  application. 


Figure  16a. 
Surface 


V- ,  Vft 
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Paraboloid  and  One  Focal 


Figure  16b.  Parapoloid  and  a  Second 
Focal  Surface 


Figure  16c.  Paraboloid  and  Composite  of  Focal  Surfaces 


-1  Surfaces  of  a  Paraboloid  Receiving  a  Plane  Wave  20°  Off-Axis 


If  it  is  desired  to  receive  the  energy  of  a  converging  wavefront  by  means  of  a 
point  source,  the  optimum  position  for  the  point  source  will,  in  general,  be  either 
on  one  of  the  two  focal  surfaces  or  between  them.  Regions  where  the  focal  surfaces 
are  close  together  or  intersect  generally  have  high  energy  density.  Degenerj.te 
focal  surfaces  are  also  regions  of  high  energy  density.  Since  every  ray  is  tangent 
to  each  of  the  two  focal  surfaces,  the  total  energy  in  the  wavefront  can  be  collected 
by  locating  properly-phased  receivers  on  one  or  the  other  or  on  parts  of  both  of  the 
focal  surfaces. 

For  purposes  of  low  side  lobes  in  a  focusing  system  there  generally  is 
a  power  taper  across  the  wavefront  whether  considered  in  the  focal  region 
or  in  the  aperture  plane.  This  taper  places  greater  importance  on  the  rays  through 
the  center  of  a  focusing  system— that  is,  rays  near  the  chief  ray— than  on  rays  near 
the  edge  of  the  system. 

In  general,  the  existence  of  two  distinct  focal  surfaces  in  the  focal  region  of  a 
focusing  system  indicates  the  presence  of  astigmatism.  For  each  different  direc¬ 
tion  of  incoming  energy,  the  separation  of  the  focal  points  and  Q2  (see  Figure  14) 
along  the  chief  ray  is  particularly  important  because  it  is  a  measure  of  the  astigma¬ 
tism  in  a  region  of  high  power  density.  In  most  one-dimensional  scanning  systems, 
the  plane  of  scan  is  a^plane  of  symmetry  of  the  system  and  also  a  principal  surface 
of  all  ray  congruences  in  the  focal  region.  The  other  principal  surfaces  for  all 
chief  rays  in  the  plane  of  scan  will  be  normal  to  the  plane  of  scan.  As  the  direction 
of  the  incoming  energy  changes  in  the  plane  of  scan,  the  two  focal  points  associated 
with  the  chief  ray  describe  two  curves.  Rays  in  the  plane  of  scan  near  each  chief 
ray  will  focus  on  one  of  these  curves,  known  as  the  T  (tangential)  curve,  while  rays 
near  the  chief  ray  but  in  a  plane  normal  to  the  plane  of  scan  will  focus  on  the  other 
curve,  known  as  the  S  (sagittal)  curve.  The  separation  of  the  S  and  T  curves 
is  thus  a  measure  of  the  astigmatism  along  the  chief  ray  for  directions  of  incoming 
energy  in  the  plane  of  scan  and  hence  is  an  indication  of  the  scanning  capability  of 
the  system.  If  the  S  and  T  curves  are  close  together,  beam  scanning  can  be  ac¬ 
complished  by  point-source  feeding  along  a  mean  curve  between  the  two.  If  multi¬ 
ple  beams  broad  in  the  plane  of  scan  but  narrow  in  the  other  dimension  are  de¬ 
sired  then  point -source  feeding  along  the  S  curve  will  produce  the  desired  result. 

If  multiple  beams  narrow  in  the  plane  of  scan  but  broad  in  the  other  dimension 
are  desired, then  point-source  feeding  along  the  T  curve  will  produce  the  desired 
result. 
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Appendix  A 

Principal  Normol  Radii  of  Curvature,  Principal  Directions,  Principal  Planes, 

and  Lines  of  Curvature 


Let  P  be  a  point  on  a  surface  W  and  let  PN  be  a  line  segment  normal  to  W 
at  P  (see  Figure  Al).  Each  normal  plane  tnrough  PN  intersects  W  in  a  curve 
C  and  each  curve  C  has  a  unique  radius  of  curvature  p  at  P.  In  general,  for 
each  point  P  there  is  also  a  unique  position  of  the  normal  plane,  such  that  the 
intersection  curve  C  =  Cj  has  a  radius  of  curvature  p  -  at  P  that  is  maximum 
(see  Figure  A2).  Similarly,  there  is  a  unique  position  of  the  normal  plane,  such  that  the 
the  intersection  curve  C  =  C2  has  a  radius  of  curvature  P  =  Rg  that  is  minimum. 

These  extreme  values  and  Rg  are  called  the  principal  normal  radii  of  curva¬ 
ture  of  W  at  P;  the  directions  of  Cj  and  Cg  at  P  are  known  is  the  principal  direc¬ 
tions  of  W  at  P;  and  the  normal  planes  whose  intersect  with  W  produce  Cj 
and  C2  are  called  the  principal  planes  of  W  at  P  .  A  ve  on  W  that  at  each 

point  has  a  direction  corresponding  to  a  principal  direction  of  W  at  tnat  point  is 
known  as  a  line  of  curvature  of  W.  It  can  be  shown  that  the  principal  planes  at  each 
point  are  mutually  orthogonal  and  hence  the  two  families  of  lines  of  curvature  form 
an  orthogonal  curvilinear  coordinate  system  on  W  . 
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_  Appendix  B 

Geometrical  Optics  Power  Flow  in  a  Source  Free,  Nonconducting, 

Isotropic,  Homogeneous  Medium 

The  ray  paths  for  this  case  are  straight  lines  (see  Section  2).  In  Figure  Bl,  P 
is  a  point  on  a  wavefront  W,  P1  is  the  corresponding  point  on  the  wavefront  W*, 


and  d  is  the  distance  from  P  to  P'  . 
The  differential  element  of  area  dA  on 
W  is  centered  at  P  and  the  correspond¬ 
ing  element  dA'  is  determined  by  the 
intersection  of  the  ray  tube  through  dA 
with  the  wavefront  W  ,  Tne  curves 
Cj  and  Cg  through  P  are  lines  of 
curvature  of  W  and  hence  are  mutually 
orthogonal  (see  Appendix  A).  The  cor¬ 
responding  curves  Cj  ai  '  C,J  in  W 
are  also  orthogonal  lines  of  curvature. 

If  the  principal  normal  radii  of  curva¬ 
ture  of  W  at  P  are  Rj  and  Rg  and 
the  centers  of  curvature  are  located  as 
shown  in  Figure  Bl,  the  principal  normal 
radii  of  curvature  of  W'  at  P'  are 
Rj  =  Rj  +  d  and  R£  =  Rg  +  d  .  Thus, 
dA  is  proportional  to  RjRg  and  dA1 
is  proportional  to  RjR£ 

If  S  is  the  power  density  at  P  and  if 
S'  is  the  power  density  at  P'  , 


Figure  Bl.  Geometrical  Optics 
Power  Flow 
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(36)  becomes 

SdA  -  S'dA' 


or 

c,  -  SdA.  .  „  R1  R2 _ 

S  dA«  '  &  <R1  +  d)(R2  +  d)  • 

The  Gaussian  curvature  K  of  W'  at  P'  is  defined  to  be 


(Bl) 


K _ 1 _ _ L . 

R{  R’  '  (Rj+  d)(R2  +  d) 

and  hence  is  a  function  of  d  .  Since  K<0)  =  ,  we  can  write 

R1  R2 

c,  -  oMdL 

b  b  K.u)  • 

This  equation  relates  the  power  density  S'  at  the  point  P'  on  a  ray  to  the  power 
density  S  at  a  reference  point  P  .  The  power  flow  so  expressed  is  equivalent  to 
the  pov;er  flow  implied  by  the  field  intensity  relation  (26). 
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magnetic  theory  in  the,  limit  as  X  -*■  0.  to  geometrical  dptiys-is  made  evident. 
The  application  of  the  "geometricaToptics  approximatiort^To  phase  analysis 
and  synthesis  is  described  and  an  example  of  synthesis  m  presented.  The 
concept  of  power  flow  in  ray  tubes  is  used  to  obtain  approximations  to  power 
distributions  in  the  antenna  aperture,  in  the  focal  reghfcn,  and  in  the  far  field. 
Ray  analysis  is  used  to  determine  those  feed  locations  in,  the  focal  region 
that  will  most  nearly  collimate  the  far -field  rays  that  If*  in  certain  desirable 
planes.  The  Theorem  of  Mai  us  is  used  to  formulate  the  <*jual  path  length  law 
and  applications  are  given.  Focal  surfaces  (or  caustics)*? elatiye  to  a  recti¬ 
linear  congruence  are  defined  and  then  used  to  present  a  .geometrical  optics 
description  of  the  focal  region.  The  equations  of  tne  fccql  surfaces  of  a  para¬ 
boloid  receiving  a  plane  wave  20*  off-axis  are  calculated  and  photographs  of 
three-dimensional  models  of  the  focal  surfaces  are  shown. 
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